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THE LEVEL DENSITIES OF RANDOM MATRIX UNITARY 
ENSEMBLES AND THEIR PERTURBATION 
INVARIABILITY 

WANG ZHENGDONG 1 AND YAN KUIHUA 1 - 2 


Abstract. Using operator methods, we generally present the level den¬ 
sities for kinds of random matrix unitary ensembles in weak sense. As 
a corollary, the limit spectral distributions of random matrices from 
Gaussian, Laguerre and Jacobi unitary ensembles are recovered. At the 
same time, we study the perturbation invariability of the level densities 
of random matrix unitary ensembles. After the weight function asso¬ 
ciated with the 1-level correlation function is appended a polynomial 
multiplicative factor, the level density is invariant in the weak sense. 


1. Introduction 


In classical quantum mechanics, the statistical properties of energy lev¬ 
els can be described by the /c-level correlation functions defined as (see 
Mehtaf]2J) 


R k 


5<X, X 2 , ■■■ , x k ) = /"■/ P n p (® 1 , ®a, ••• , x n )dx k+1 ■■■dx„ 


where P ri/3 (x 1 , x 2 , — • , x n ) = c n/3 -exp(— f3H) is the joint probability density 
function of n eigenvalues x x , x 2 , ■ ■ ■ , x n of a n x n random matrix, c n/3 is the 
normalized constant. H is the Hamiltonian of the logarithmical interacting 
n particles system on a straight line, which is given by constraining one-body 
potential and logarithmic repulsive two-body potential, i.e. 


n 

H = ^2 v ( Xi )- X iogl^-x-l, arGK. 
2—1 


In general, V(x ) is called potential function and (3 Dyson’s index. f3 = 1, 2 
and 4 are corresponding to the orthogonal, unitary and symplectic ensembles 
respectively. When k = 1, R^(x) can be explained as the distribution 
density of energy levels which can be found near by x. The level density 
denoted by o^ (x), which is a global quantity, is defined by the limit of 
the 1-level correlation function R l n p(x). Then how to determine the level 
density? It can be traced back to Wigner’s pioneering work |2()l 121 j . The 
results of early work are reviewed in da ce]. Recently, there are many 

l 



authors to concentrate on this problem (See SpohmTYj, Bai and Yinj^j, 
Nagao and WadatipHj. Haagerup and Thorbjprnsen(Hj. GirkojB], Kiessling 
and SpohniTU, Duenez LedouxjTTj , etc.). 

Notice that given different or special one-body potential V(x), it will 
exhibited kinds of images for us. As a matter of fact, in the case of classical 

x 2 

Gaussian ensembles, V(x) = —, x £ R. The level density is the famous 
“semicircle law” first derived by Wiener pi PIT . i.e. 


cr (*) 


\'j2n — x 2 x 2 < 2n 
0 x 2 > 2 n. 


In addition, in the case of Laguerre ensembles, the level density can be 
evaluated by a physical argument (See BronkJ^), i.e. 

( —7= \/2n — x 0 < x < 2n 
oJx) = < 7r v i 

pK 1 \ 0 x > 2n. 


In the case of Jacobi ensembles, the level density can also be evaluated 
by a physical argument (See Leff|12jh i.e. 


a^x) 


n 

7T\/1— X 2 

0 


— 1 < x < 1 
otherwise. 


In the case of unitary ensembles, i?,) 2 (x) can be expressed to a concise 
formula (see uni, d) which is closely correlated to classical orthogonal 
polynomials, i.e. 


71—1 

(!) R L( X ) = J2p 2 m( X )™( X )’ 

771=0 

where w(x) = cexp(— 2V(x)), c is the normalized constant and p m (x) be the 
m-order normalized orthogonal polynomials associated with the normalized 

weight function w(x), i.e. J p m (x)p n (x) ■ w(x)dx = 5 mn . 

In U], Haagerup and Thorbjprnsen only studied the Gaussian unitary 
ensemble(GUE) which was denoted by SGRAl(n,o 2 ) there. Using the fol¬ 
lowing property of Hermite polynomials 

k 

H k (x + a) = J2 C k( 2a ) k ~ j H j (x), aeM, 

3=0 
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the authors directly obtained an equality for the complex Laplace transform 

of - R L( x ^ le ■ 

n 

/* 1 /* 1 ^ ^ 

exp(sx) i^—R l n2 dx = J^exp(sx) £ R l( x )e~ x2 ^j dx 

2 2 

= exp(^-)<L(l - n, 2; -y), s£ C, 

where H k (x ) is the k-order normalized Hermite polynomial, ^(a, c; x) is the 
confluent hyper-geonretric function(CHGF) with parameters a and c. Then 
they gave a short proof of Wigner’s semicircle law. 

As we know, the CHGFs are complicated series expansions. In nu, 
Ledoux pushed forward the investigation by Haagerup and Thorbjprnsen 
and only concentrated on the differential aspects of CHGFs. The author con¬ 
structed an abstract framework of Markov diffusion generators, and in which 
derived the basic differential equations on Laplace transforms of p‘^ n (x)w(x). 

Using the recurrence formula of Hermite polynomials and the uniform 
integrability of random variable sequence, by the obtained differential equa¬ 
tion, the author showed that 

m— 1 

f(lTT=)— £ PfcOVO) dx = E(f(VXY)\, for all / G C b ( M) 

which determines the level density of GUE in the weak sense, where X and Y 
are two independent random variables with the uniform distribution on [ 0 , 1 ] 
and the arcsine distribution on (—1, +1) respectively. By the analogous 
technique, the author also obtained the level densities of Laguerre and Jacobi 
unitary ensembles respectively. 

In this paper, we will generally deal with the problem for kinds of unitary 
ensembles (i.e. (3 = 2) by operator method. It is no confusion to omit the 
subscript 2 in the below text. Moreover, as we will see, this method can 
effectively be used to study the perturbation invariability of level density. 

It is well known that the normalized orthogonal polynomials p n (x) satisfy 
the following recursion formula (see section 2 for details) 

xp n (x) = a n p n+1 (x) + p n p n (x) + 7nPn-i(*)> n = 1, 2, 3, ■ ■ ■ . 

In the paper, we assume that a n and f3 n satisfy the following exponential 
growth conditions: 

( 2 ) «»=^(l+0< P n = Cn\l + C n )+ Vn, 

where £/ 0 >C;i 0 > 0 <t<l are constants and lim £ n = lim ( = 

n —>oo n—>oo 

lim rj n = 0 . 
n—>oo 
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Note that the classical Hermite, Laguerre and Jacobi polynomials all sat¬ 
isfy exponential growth conditions (see section 2). 


We define the “ascending”, “equilibrating” and “descending” operators as 
follow, 


A +Pn{ X ) = «nPn+l(®), A 0 Pn( X ) = PnPn ( X ) > A -Pn ( X ) = TnPn -1 ( X ) ■ 


In section 3, we use these operators to obtain the moment inequality of 

the probability density cr n (x) = R\{x-sJ D n ) and to show that the limit 

of k-th moment of cr n (x) exists (see theorem 13.211 . where D n is the 2nd 

moment of — R 1 (x). Further more, by considering the limit of characteristic 

n " 

functions of a n (x), we obtain the main result in section 3, i.e. if there exists 
a probability density cr(x) with k-th moments M^ = lim M^ k \ then (see 

n— >oo 71 

theorem 

n —i oo, 


where w means weak convergence. 


In section 4, the perturbation invariability of level density is studied. 
We will consider another weight function w(x) by appending a polynomial 
multiplicative factor p 2 (x) to w(x), which uniquely determines a family of 
orthonormal polynomials {p m (x)}. Accordingly, 1-level correlation function 

n— 1 

is Rl(x) = E p 2 .{x)p 2 {x)u j(x). We will show that the probability density 

3=0 

&n( x ) = R\(x\JD n ) is still weakly convergent to a(x). 


2. Orthogonal Polynomials 


In this section, we briefly introduce the orthogonal polynomials. We 
would like to lead the readers to refer 0EEUIE1. Let us start by quot¬ 
ing a classical result as follow. 


Theorem 2.1. The following relation holds for the normalized orthogonal 
polynomials p n (x) with normalized weight function m(x): 


xp n {x) = a n p n+1 {x ) + P n p n (x) + 7„P„-i(»), n = 1, 2, 3, • • • 


where a n , (5 n and r y n are constants which can be expressed in terms of the co¬ 
efficients a n and b n of the highest terms inp n {x) = a n x n -|-6 ri a; n ” 1 • • • , a n ^ 
0, i.e. 


a„ = 


n+1 


Pn = — ~—, n = 0, 1, 


n+1 


7o = 0, n = 1, 2, • • • • 


In = «„-l, 
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If we regard the multiplication by x as an operator A x , clearly the operator 
A x plays “ascending”, “equilibrating” and “descending” roles when it acts on 
the n-order orthogonal polynomial p n (x), i.e. 

A x = A + + A 0 + A_, 

where A + , A 0 and A_ are called “ascending”, “equilibrating” and “descend¬ 
ing” operators respectively and defined by 

A + p n (x) = a n p n+ i(x), n = 0, 1, - ■ - , 


(3) A oPn (x) = P n p n (x), n = 0, 1, 

A_p n (x) ='Y n p n _ 1 (x), A_p 0 (x) = 0, n = 1, 2, ■ ■ ■ . 

In the below text, we always assume that a n and (3 n satisfies the expo¬ 
nential growth conditions ©. Let us consider several classical examples. 

1) Her mite polynomials 
The weight function 

w(x) = exp(— x 2 ), x € R 

and 


Hjx) = 


1 


V2 n n\ 


H n {x), 


where 


d n 


H n( x ) = (-l) n exp(x ) — (exp(-x )) 


is the standard n-th Hermite polynomial which satisfies the orthogonal re¬ 
lation 

TT~ 2 [ H k (x)H j (x) ■ w(x) dx = 2 k k\S kj . 

J R 

H n (x) satisfies the following recurrence relation 


xH n {x) = \/-y-^n+i(*) + J^H n -iix), n = 1,2, 3, 


H 0 (x) = 1, H 1 (x) = V 2 


X. 


/ Jl I I 

P n = 0, i.e. f = -J=, C = o, t = i ri = C n = rj n = 0. 

2) Laguerre polynomials 
The weight function 

w^ a \x) = x a e~ x , x 6 (0, oo) 


U a \x) = a/ tv ^ + 1) , -L (a) (a), 

" |r(Ha + i) " 


and 








where 


ry — apX pin 

L&(x) = , r ( X n+a e~ x ) 

n n\ dx 11 

is the standard n-th Laguerre polynomial which satisfies the orthogonal 
relation 

jf li‘ ) (# ) W -^e-'dr = r( r ( t + t) 1> ^- 

(x) satisfies the following recurrence relation 
xL\f>{x) = - v /(n + l)(n + l + a)L^ 1 (x)+(2n+a+l)L^ (x)-^n(n + a)L(°_\(a;), 

1 


L(“ } (x) = 


__=, tW(.)= ~* + ° +1 , „ =1,2,3,.... 

vT(a+l) v ' r (“+ 2 ) 


Thus a n = - y/(n + l)(n + 1 + a), /?„ = 2n+a+l, i.e. £ = —1, C = 2, i = 1, 
c _ , 2 + a 1 + ao /2 iA_^ + a _ n 

C - (H --f- C„ - ’ dn - 0. 


2 n 


3) Jacobi polynomials 
The weight function 


and 

where 

and 


w^ b \x) = {l-x) a (l+x)\ z€(-l, 1) 

J(“’ &) (z) = ( c (“’ b ))-§j(“’ & )(x), 




( a , M _ 2 a+b+1 T(fc + a + 1)T(A; + 6+1) 

1 T(fc + 1 

(-l) n d" 


Cfc 2fc + a + 6+ir(/c + l)T(fe + a + 6+l) 


1 


■((1 -x) n+a (l + x) n+b ) 


(1 - x) a {l + ,x) 6 2 n n! dx r 

is the standard n-th Jacobi polynomial which satisfies the orthogonal rela¬ 
tion 

£ Ji a ’ b \x)j^ b \x) • (1 - x)°(l + x) b dx = c^S kj . 
j( a,b )(x) satisfies the following recurrence relation 


xJ^ b \x) = 


(n + l)(n + a + l)(n + b + l)(n + a + b + 1) 


2n H - CL H - b H - 2 


(2n H - cl b -|- 3)(2n -|- a -(- 6 H- 1) 


J ( “’ b) (.x) 

n+1 ' / 


a 2 — 6 2 


+ 


(277/ CL ~\~ 6) (277/ ft -|- b 2) 
2 


j(a,fe)( a 

n 


/ 77 /( 77 / + a) ( 77 / + 6) ( 77 / + a + b) 

2n cl by (2n -h ci -|- 6 H - 1)(277/ -|- a H- 6 — 1) 


J (a ’ b) ( a 

n—1 v 


ThUS f = 2’ C = °’ * = °- C “ d h ” e th6ir reSpeCt ™ eXpreSSi0nS ' 



In particular, if a = b = 0, then it is the Legendre polynomial. The 
recurrence relation is 


*j(°>°)(x) = 


n + 1 


n 


_=j(°’°)(x) + 

^(2n + 3)(2n + 1) n+1 ^(277 + l)(2n - 1) 




3. The Level Densities 


3.1. Main Results. In this section, we consider the limit of the global 1- 

n— 1 


level correlation functions /?,’ (x) = p 2 m {x)w{x) for unitary ensembles, 


m =0 


where p m (x) is the m-th normalized orthogonal polynomial associated with 
the normalized weight function w(x) = cexp(— 2V(x)), c is the normalized 
constant. 


Specially, for the Gaussian unitary ensemble (GUE), V{x) 


x 2 


m(x) = 7T 2 w(x), 


Pm(x) = x G M, 


where H m (x) is the m-th normalized orthogonal Hermite polynomial asso¬ 
ciated with the weight function w(x) = exp(—x 2 ). 


For Laguerre unitary ensemble (LAUE), V (x) 


—a log x + x 0 <x<oo 
oo otherwise ’ 


w(x) = cw { - 2a \2 x), p m [x) = v / 27 cL^ a) ( 2 x), x € (0, oo), 


where (x) is the m-th normalized orthogonal Laguerre polynomial asso¬ 
ciated with the weight function w ( ' a \x) = x a e~ x . 

For Jacobi unitary ensemble (JUE), V (x) = 


— log((l — x) a (l + x) b ) —1 < x < 1 
oo otherwise 


zu(x) = cw {2a ’ 2b \x), p m (x) = \ZTfcJ^ a ' 2b) (x), x G (-1, 1), 

where J^’ b \x) is the m-th normalized orthogonal Hermite polynomial as¬ 
sociated with the weight function w^ a ’ ^ (x) = (1 — x) a (l + x) b . If a = b = 0, 
it is the Legendre unitary ensemble (LEUE). 



Now it is no less of generality to consider the whole real axis K. Denoted 

by D the 2nd moment of the probability density — R}(x). Then 

" n " 


n— 1 


If 1 1 

D n = - x 2 Rl(x)dx = - ^2(x Pj , x Pj ) 

n n j =0 

- n— 1 

= ;Etf+tf+7?) 


3=0 


But 7 . = a._ 1 , hence 


n —2 


(4) 


A* “ “ + + “(“Li + ^n-l)‘ 


i=o 


Let 

^n(^) = 

Then the k-th moment A/^d of probability density cr n (x) is 

Af ( fc ) = f x k a n (x) dx = - 1 , [ x k P 2 (x)w(x) dx 

Jr ” n-(D n )Wj^J R y ’ 


71—1 
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^((^ + +^o+^-)S>P J 


n ■ ( D n) k/2 —“q ' ' LH*,) 


Let A* be an operator set composed of those terms in the expansion of 
{A + + A 0 + A _) k , in which the operators A + and A_ exactly appear i times. 
And note that for all T 0 I I A* (Tp., P .) , =0, thus we obtain the 

i 

following lemma 

Lemma 3.1. 


n— 1 L 2 J 


(5) 


Al (fc) = 


n 




EE E (iK. 


j=0 *=o TeA* 


r 2 (ro) 



I 2-I 

Moreover, set r^k) = C* -2i , r a (fc) = (2£ 2 + £ 2 ) fc/2 and 

i =0 

(fc) (2t + l) fc / 2 

r 2 (fc) (fct+l) 

then we have the following limit theorem for the k-th moment M^ k \ 
Theorem 3.2. Under the exponential growth conditions, 

(6) lim = M^ k \ for any k G Z + . 

n—>o o n 

The proof is based on elaborate estimations of M^. Here we firstly 
consider the classical cases instead of being anxious to verify this theorem. 


1 1 C m 

Example. For GUE, £ = —, £ = 0, t = so ih/( 2m ) = -2s-, M^ 2m+1) = 

V2 2 m + 1 

0. Moreover, M^ is exactly the k-th moment of density function 

, 2 

x 2 > 4. 


^gO) = <! 2u 


1 \/4 — x 2 x 2 < 4 


For LAUE, £ = -1, £ = 2, t = 1, so M (fc) = 




2 fc /2(fc + i)' 


Moreover, 


V2V2- 

tt\/2x 


X 


0 


Af( fe ) is exactly the k-th moment of density function 

0 < x < 2\/2 
x > 2 1 / 2 . 

1 pm 

For JUE, £ = C = 0, t = 0, so M (2m) = -2 e-, M (2m+1) = 0. Moreover, 
MW is accordingly the k-th moment of density function 

<Tj(x) = — , „ , — 1 < x < 1 . 


7T 


\/l — X 2 


Now suppose that there exists a probability density &(x) with k-th mo¬ 
ments M^ k \ Then we have 

Theorem 3.3. Under the exponential growth conditions, 

(7) a n (x) a(x), n-^ 00 . 

where w means weak convergence. In general, a(x) is called a level density. 


Corollary 3.4. The level densities of Gauss, Laguerre and Jacobi unitary 
ensembles are <j g (x), 0 L (x) and (Jj{x) respectively in weak sense. 

By the above example and theorem IQ it is obvious. 


9 












3.2. Proof of main results. Now we begin to verify the main results. First 

of all, let us consider the 2nd moment D of —R 1 (x). If C > 0, set 

n n 


( 8 ) 


u. = max 

3 j—k<m<j+k 


{iCJ, ICJ + > t }> 


(m t 


then lim u. = 0. It is no less of generality to assume u. < 1. Thus by the 


3^o o 


exponential growth conditions a and equality a, 

2e+c 


(9) 


n 


^/‘(l - < 0„ < ?L±±. + 


3=0 


3=0 


If £ = 0, one can easily choose another proper infinitesimal sequence u j such 
that the above inequality is still valid. By a, it is easy to see that 

( 10 ) 


3 = 0 ( 1 ). 


n 


Proof of theorem \S.^[ Firstly, we come to verify a result about series limit. 

oo 

That is for any sequence x n and y n . if lim x n = 0 and y n > 0, > y n = oo, 

n—>oo J 

then 


n=0 


( 11 ) 


lim 


E»: 

3=0 


> x i 


3=0 


= 0. 


Indeed, Me > 0, 3N E N, as n > N, \x n \ < |. Thus 

n— 1 N n 

Y yj\ x >\ Y y : , 


3=0 


< 


3=0 


— n 


Y y> Y Y y. 


£ j=N £ £ _ 

+ 2 ' ~JL “ 2 + 2 _ c ' 


3=0 3=0 j =0 

Now let us set two cases to verify the theorem. 

Case 1. C > 0. 

By exponential growth conditions a and the definition 0 of u j , we see 
that 
(! 2 ) 

,max {|aj} < |£|(j+£O t (l+u ), . max, {/3 m } < +kf (l+u ). 


j-k<m<j-\-k 


j—k<m<j+k 
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Thus for j > k and T G A*, 

{Tp >’ p >> s C- k ?s + d |a ~ |> ) 2 ‘C-»§Si + d /3 " > ) 

< (|£IW + *)‘(1 + »,)) 2 '(ai + *)‘(1 + <*,-)) 

And then 

(13) (Tp v p^<e i C, k - 2 \j + k) kt { l + Ui ) k . 


k—2i 


k—2i 


Therefore by lemma mi and inequalities (THT1) .(E)1)- (I1()I) . 


n —1 


(14) M (k) < ^- J — 


^ j + k ) kt (l+‘ 


r 2 (k) ^ n “ 2 

n 


n—z 

71 L ^ J 


k/2 


+ 0 (Tfcm)> forfcGZ+. 


n 


j =o 


On the other hand, also by exponential growth conditions © and the 
definition © of u ., 

(15) 

. min {|aj} > ItKj-kf min {/3 m } > (/{j-kfil-u^. 

j-k<m<j+k j—k<m<j+k 


Then for j > k and T G A*, 

J k 7 

< Tf V p 4 a C- t |Sii + d |a " l >) J 'C-^<j + d /J " } ) 

> («i-fc)‘(i-'«,)) 2 ‘(c<i- t )‘( 1 -“ J )) 

Therefore, 

(16) (TV, p.) > e 2i C fc - 2 *0' - k) kt ( 1 - u,) k . 


k—2i 


Thus by lemma mi and inequalities © - (11 Of) and m, 


71—1 


(17) M {k) > ^ 

" _ r 2 (/c) ’ n ~ 1 


1 x 

n ' J 


1 71—1 

-5^j 2t ( i + 

n r—' 


u. 


k/2 


+ 0 (zm+ i)> for n > 


n 


3=0 


li 








Now let us come to consider the limit of the dexter series in inequalities 
<d and m- By equality CD, 


^ n—1 

-£o+*)“(!+■ 


(18) lim 


1=0 


1 


n 


n—1 


£(j + *0“ 


n —>oo , j n—2 


= lim 


l=o 


(2f + l) fc/2 




sk/2 n ^°° ( 1 oA fc/2 (fct + 1) 


1=0 


- y^ j 2< 


and 


(19) lim 


n—1 


1=0 


n—1 


n 


£(i- *)“(! -«j)' 


!=*0 


n 




kt 


n —>oo , 2 n— 1 


= lim 


l=o 


(2i + l) fc/2 


~X] j ' 2, o+“j) 

n J 


k/2 


n >oo , i n—2 fc/2 (Atf + 1) 

3 


n—z 

_ „'2i 

n 


1=0 1=0 
Then by inequalities CD and CD and equalities d-d, we have 


( 20 ) 


lim M< fc) = ■ ( 2t + 1 ^ k/ ~ = M (k) if C > 0. 

n—nx n r 2 (k) (kt + 1) 


is 


Case 2. ( = 0. 

Set v. = max j|£l, Ini}. It is obvious that lim v. = 0. So it i 

3 j—k<m<j+k'' 1 j —>oo 3 

no less of generality to assume that v. < 1. Thus by exponential growth 
conditions Q, for 0 < j — k < m < j + k, 

(21) l£l(j-&)*(!-^) < l«J < \£\ti+ k) t (l + v j ), \P m \<v r 
Of course, v. may be identity to zero. So we come to discuss the estimation 
of MW in terms of the parity of k. 

(i) If k is odd, then for j > k and T £ A*, 

(TP, , Pj ) < (fO' + AO^l + w,-)) 


2 i 


,k—2i 


<f i v k - 2i (j + k) Kt (l + «.)". 


kt ( 


Let v j = yy 2l , thus by inequality (EHl- HUl) and lemma mi 


i=0 


( 22 ) 


M 


< 


1 


1 


n 


n—1 


Eo +fc )“( i + 


v) k v- 

J 2 ^3 


1=0 


r 2 (fc) ’ n “ 2 

Vn 


.j n—z 

lET(i-»T) 


fc /2 


+ °(zfcm)’ forfceZ+ - 


n 


l=o 
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Note that lim v. = lim v. = 0. Then by equality m and CHI, 

J—>00 3 i-OO 3 


lim 

n—>oo 


^ n —1 

-£cj + *o“( 1+ >T®. 

3=0 


n —2 


a- ^ 


\ fc /2 


j=0 


= 0. 


Thus by inequalities 11221 ) . 

(23) lim = 0 = (here rAk) = 0). 

n^oo 71 


(ii) If k is even, then by inequality CH, for j > k and T £ A*, i A 
(24) 

and for T G A k / 2 , 


(T Pj , Pj ) <ev k - 2i {j+k) kt {i+ Vj y 


(25) eu - fc) fct (l - U,) fc < (Tp., p.) < eu + kT(l + v 3 y. 


Set 


and 


71—1 


-j lb X 

^ = „.m w EEEfa.-' 

n ^ j=o#kTeA* k 


P, 


71—1 


"»‘ , = 

1 j=° reA p 


p,-;> 


then by equality ©, 


M (fc) = M (fc) + 


n n 


It is completely analogous with the above discussion in the case of odd k, 
we see that by inequalities © and (121) and equality CD, 


(26) 


lim M (fc) = 0. 


But with respect to M^ k \ by inequalities © and (121 . we have 


1 


71—1 


(27) Ml* 1 < ii ■ J± -- + 0(-kt T ) 


r 2 (fc) - "- 2 
Vn 


n zJ 


fc /2 


n 


3=0 
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and 


r^k/ 2>./c Ti . 

(28) M (fc) > ^ ' 


n—1 

n ^^ J 


r 2 (k) 


1 n—2 

-Ej 2t ( 1+u J -) : 

n r — 1 ' 


j=o 


k/2 


+ o{- 


1 


n 


kt+l' 


Then by equalities m-m and inequalities dZZD - dlHD , 

a(t) = ++ . (2t+i) t/2 
n r 2 (k) (kt + l) 

Therefore by equalities m and El, 


(29) lim 


(30) 


lim M w = C T+ . = nh). < a + 1 > t/2 = M <*> 

n—>oo » r 2 (fc) (fci + l) r 2 (lfe) (fci + 1) 


Combining equalities ei, (Ei and El, we complete the proof of the 
theorem cm 

Next we come to verify the theorem 13.31 In the first place, we have the 
following rough estimation of k-th moment M^. 

Proposition 3.5. Under the exponential growth conditions for any 
e > 0 there exist an integer K 0 (s) independent on k and n such that 

M (fc) 


(31) 


<d k (ek) kt , for k,n > K 0 (e), 


where d 0 is a constant independent on k, n and e. 


Proof. By the exponential growth conditions ©> there is a constant d 1 
independent on k and n, such that 

Kl <dX, 

Then for any T € A*, 


(32) 


(T Pj , Pj ) <d k (j + k ) 


kt 


And by inequality ©, there is a constant d 2 > 0 independent on k and n, 
such that 


D n > d 2 n 2t . 


(33) 

Thus by lemma ITT! and inequalities (I32l) - (l33l) . 


n —1 


M (fc) 


< 


E c i c i-i £ d "+ k)kt d i E 

i =0 _ 3=0 

(d+ /2 


k—i 


n 


fct+1 


< 


2=0 


(d+ 2 


n + k\ kt 


n 
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But ^ CjCti < 3 fc . Set d 0 


i =0 


3ri, 

W 


, then 




M< fc > + 

" - °V nJ 


Note that for all e > 0, set K 0 (s) = [|] + 2, then when k,n > K 0 (e), 
k £ £ 

1 H— < —k + —k < ek. Consequently, 


M (fe ) 


<d k (ek) kt , for k,n> K 0 (e). 


Thus we complete the proof of proposition 13.51 

Proof of theorem IT 51 Let f n (9) and f(9) be the characteristic functions 
of cr n (x) and cr(x) respectively. It is sufficient to show that 

lim f n (6) = m. 


Note that 

00 l«| k 


°° \n\ k (• °° | n 1 2 m 00 

/ \x\ k a n (x)dx = Y P—M^ + Y 

^ k\ J 1 1 nK J ^ (2m)! " 

k =0 m =0 v ' m =0 


00 j^|2m+l 


m =0 

v \a\2m 

< y P—M 2m) + Y 

z -"' (2m)! 71 z -"' 


(2m. + 1)! 




|2m+l\ 


m=0 

00 |^|2m+l 


Tf ( 2m + !) 


-(1 + Mp m+2) ). 


Then by the rough estimation (1311) of M ^ in proposition 13.51 for any 9 we 
can choose e > 0 such that 


00 \e\ k 

~k£ 


E 

k =0 


\x\ k a n ( x ) dx < oo, for n > K 0 (e). 


Thus we have 


/..Mi = J e' ex <r„(x) dx = J jf ^ff-ojx) dx 


E 


( i9x ) 


fc =0 

OO 


fc! ^ (*) dx = Y ^r M P > for n ( £ ) ■ 

fc =0 


_ (iQ\h 

Analogously, it can be verified that f{9) = Y, M^ k \ 

k =0 

Note that 0 < t < 1. So for any 0 we can choose e > 0 such that 

g |gd 0 l fc (efc) H < ^ 


k=K~ (e) 


fc! 
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Therefore by theorem IM.21 inequality (Kill) and Lebesgue control convergent 
theorem for series, 


lim 

n—>oo “ 


V = V lim = V 

' fcl " 4— i k\ n-> oo " ^ 


(*0) A 

“fcT 




fc =0 


k=0 


k =0 


Thus we obtain 

lim / n (0) = /(0). 

n->oo 

It completes the proof of theorem 13.31 

4. Perturbation Invariability 


Now Let p(x) be a fixed l-order polynomial, denoted by p j (x) the j- 
th normalized orthogonal polynomials associated with the weight function 
w(x) = p 2 (x)m(x). By equality 0 , the 1-level correlation function is 

n —1 

#(®) = ^2p 2 j( x )p 2 ( x ) w ( x )- 

3=0 

Let 



Set 

H n = span{p 0 (x ), p x (x), ■ ■ • , P^x)}, 

then H n is a subspace of L 2 (M, zu(x)dx) with n dimensions. It is obvious 
that {p 0 (x)p(x), ■■■ , p n _ l _ 1 (x)p(x)} is a family of normalized orthogonal 
vectors in H n . We can extend this set of vectors, such that it makes up of 
a normalized orthogonal base of H n . Denoted it by 

{e { ^\x), ••• , e[5(x), p 0 (x)p(x), ••• , (x)p(x)}. 

Let P n be the projective operator from L 2 (R, w(x)dx) to H n . we con¬ 
struct an operator from H n to itself as follows, 

(34) T {k) = P n oA k , 

where A x is multiplication by x. 

Now we consider the fc-moment M^ and M^ of probability density 
ar (x) and <r n (x) respectively. 


Proposition 4.1. 

(35) 


M (k) 


Tr(TW) 

n-(D>/2' 
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Proof. 


M (fc) = / x V n (x)dx = 


n —1 


n-{D n ) k / 2 j^ Q 




L 2 (U7) 


1 


n —1 




L 2 (td) 


1 


n—1 

E 


n ■ ( D n ) k/2 ^ n ^ 


Tr(T < fc) ) 

' ; ’ P ^, 2 ( ro) “n-(D„>/2- 


Lemma 4.2. 

(36) 


ll^-/ILa( w ) < 3JVJI/II, forallfen n , 


where N n = sup {|aj, |//|}. 
0<i<n—1 


n— 1 


Proof. Put / = ^^c i p i , then 


i=0 


n—1 


ll^/ll 2 = II Will 2 = HlZ c i( a ^+i + +7iPi-i) 

i=0 * 

= ll5Z( c i-i a i-i + C *A + c i+i7i+i)Pil| 2 

i 

= X]l c i-I a i-1 + c iA +c i+ iaj 2 

i 

< 3N 2 V (c 2 , + c 2 + c 2 ) < 9iV 2 

— n / v v i— 1 % i+l ; — n 

1 

Therefore, ||AJ|| < 3iV„ 

Corollary 4.3. 


fc-i 


(37) 


l^/ll<3 fe (n^ + ,)ll/ll, /or all f €E 7i n 


3=0 


Proof. By the above proposition, it is obvious. 

Proposition 4.4. 

(38) M (fc) = + I, 

where 


^ n —1 /—I 

^mit^E -EW).^ , 

v«/ j=n—l r ( ro ) j=o r^O) 
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Proof. 


Mf } = [ x k a n (x)dx 
JR 






n—l—1 


n • (D n ) k / 2 V 


E 



+ 


n 


1 n—1 

, [D \k/2 ( E PjP 

V j=n—l 


l 2 0) 


Tr(TW) 

n-(E) fc/2 


+ / • 



Then by proposition 14.11 we obtain the conclusion. 

Theorem 4.5. Under the exponential growth conditions, 
(39) lim M^ = lim M ^, for any k € Z + . 


Proof. By the exponential growth conditions and corollary 14.31 we see that 


there exist constants C 1 , C 2 > 0 such that 


(40) 


and 


A k Mv), p,p) I < IM^pII < s‘Jv‘ +1+l _ a < c,»“ 

' r 2 /_ 


L 2 (U7) 


(41) |(Tf fc )(eW), ef)) | < ||PJ| • ||A* e W|| < 3 k N k +k _ i < C 2 


n 


kt 


i 2 (ro) 


Thus by proposition 14.41 as n » 1, 

n—1 


\i\ < 


i 


i-i 




E C 1 n kt + '£ i C a 


n 


kt 


3=0 


{C 1 + C 2 )l n kt 1 

n (D„) k / 2 n 
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Whereupon we obtain 


lim = lim M®. 

n—>■ oo n n—>oo n 

Remark. In the case of finite interval (—s, s) C R, note that the operator A x 
is bounded, so the inequalities m and m in the above proof are naturally 
valid. 

Finally we have the following theorem which tells us that after the weight 
function m(x) is appended a polynomial multiplicative factor, the limit be¬ 
havior of the normalized 1-level correlation function is unaffected in the weak 
sense. 

Theorem 4.6. Under the exponential growth conditions, 

n —> oo, 

where w means weak convergence. 

Proof. It is completely analogous to the proof of theorem Id.Ml here we omit 
it. 
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